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A path integral formulation is developed for the dynamic Casimir effect. It allows us to study 
small deformations in space and time of the perfectly reflecting (conducting) boundaries of a cavity. 
The mechanical response of the intervening vacuum is calculated to linear order in the frequency- 
wavevector plane, using which a plethora of interesting phenomena can be studied. For a single 
corrugated plate we find a correction to mass at low frequencies, and an effective shear viscosity 
at high frequencies that are both anisotropic. The anisotropy is set by the wavevector of the 
corrugation. For two plates, the mass renormalization is modified by a function of the ratio between 
the separation of the plates and the wave-length of corrugations. The dissipation rate is not modified 
for frequencies below the lowest optical mode of the cavity, and there is a resonant dissipation for 
all frequencies greater than that. In this regime, a divergence in the response function implies 
that such high frequency deformation modes of the cavity can not be excited by any macroscopic 
external forces. This phenomenon is intimately related to resonant particle creation. For particular 
examples of two corrugated plates that are stationary, or moving uniformly in the lateral directions, 
Josephson-like eff'ects are observed. For capillary waves on the surface of mercury a renormalization 
to surface tension, and sound velocity is obtained. 



I. INTRODUCTION AND SUMMARY 

The Casimir effect is a surprising phenomenon in which quantum fluctuations of the electromagnetic field in 
the vacuum, subject to the boundary conditions imposed by two conducting plates of area A and separation H, lead 
to an attractive force between them, given by 

F^te.tic{H)^-—xhcj^. (1) 

This force, which provides a direct Hnk between quantum field theory and the macroscopic world, has been measured 
experimentally (using a torsion pendulum) between a gold plate, and a gold plated sphere, to the accuracy of within 
%5 . Related fluctuation-induced phenomena occur in various areas of physics ranging from cosmology to statistical 
mechanics of phase transitions [^|J^] . 

Since the pioneering work of Casimir, various modifications to the standard boundary conditions have been studied. 
One generalization, relevant to experimental measurements of the force, is to introduce deformations (roughness) of 
the surfaces. For example, in Ref. |^] a multiple scattering approach is used to compute the interactions for arbitrary 
geometry in a perturbation series in the curvature. In another phenomenological approach, introduced in Ref. 
small deviations from plane parallel geometry are treated by using an additive summation of van der Waals-type 
pair-wise potentials. The former method is difficult to implement, while the latter is not always reliable due to the 
additivity assumption Q. 

New phenomena emerge for moving boundaries in the generalization to the dynamic Casimir effect [p^|-|l9|. The 
creation of photons by moving mirrors was first obtained by Moore pO| ] for a 1 dimensional cavity. Fulling and Davis 
found a corresponding force proportional to the third time derivative of the mirror displacement. Their methods 
specifically make use of the conformal symmetries of the 1+1 dimensional space time, and hence are not applicable 
for the more realistic case of 1-1-3 dimensions [ po| . 

Another approach to the dynamic Casimir effect starts with the fluctuations in the radiation pressure on a plate 
p3| , p^ . The fluctuation-dissipation theorem is then used to obtain the mechanical response function, whose imaginary 
part gives the dissipation in response to motion. This method does not suffer from the causality problems of earlier 
approaches |^0|, and can be applied to all dimensions. For example, the force in 1+3 dimensional space-time depends 
on the fifth time derivative of the displacement. Several other approaches have focused on the emission of photons 
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by a vibrating cavity [[15|-|17[. This is generally too small to be experimentally detectable |£8 22|; the most promising 
set-up is the resonant production of photons when the mirrors vibrate at the optical resonance frequency of the cavity 
psf . More recently, the radiation due to vacuum fluctuations of a collapsing bubble has been proposed as a possible 
explanation for the intriguing phenomenon of sonohiminescense . A good review of the topic including more 

extensive references can be found in Rcf. psf . 

A number of authors have further discussed to notion of frictional forces: Using conformal methods in 1+1 dimen- 
sions, Ref. I^J] finds a friction term 

FmctioniH) ^ a Fst^tUH) (^^^ , (2) 

for slowly moving boundaries, where a is a numerical constant that only depends on dimensionality. The additional 
factor of (v/c)^ would make detection of this force yet more difficult. There are a few attempts to calculate forces (in 
higher dimensions) for walls that move laterally, i.e. parallel to each other [^-^: It is found that boundaries that are 
not ideal conductors, experience a friction as if the plates are moving in a viscous fluid. The friction has a complicated 
dependence on the frequency dependent resistivity of the plates, and vanishes in the cases of ideal (nondissipating) 
conductors or dielectrics. The "dissipation" mechanism for this "friction" is by inducing eddy currents in the nonideal 
conductors, and thus distinct from the Casimir effect. 

In this paper we present a path integral formulation, applicable to all dimensions, for the problem of perfectly 
reflecting mirrors that undergo small dynamic deformations. Although the original Casimir problem with parallel 
plates can be tackled with such a path integral approach (see Appendix it is much easier to solve by the standard 
operator method, followed by expansions in modes appropriate to the specified geometry and boundary conditions. 
However, if one wishes to study the effect of deformations of either static (fixed roughness) or dynamic nature, the 
standard method is hardly tractable p|,p8t. In contrast, the path integral method is much better suited to handling 
deformations. The only limitation of the method is that it provides the result as a perturbative expansion in powers 
of the ratio between the deformation field and the average separation of the plates. In the present context, a great 
advantage of this method is that it does not make any distinction between "time" and "space" coordinates provided 
that appropriate Wick rotations are performed. Consequently, one can consider time dependent boundary conditions 
with no additional complications. 

Using this procedure, we calculate the mechanical response function in the frequency-wavevector domain. It is 
defined as the ratio between the induced force and the deformation field, in the linear regime. From the response 
function we extract a plethora of interesting results, some of which we list here for the specific example of lateral 
vibrations of uniaxially corrugated plates: 

• (1) A single plate with corrugations of wavenumber k, vibrating at frequencies lo <^ ck, obtains anisotropic 
corrections to its mass. Thus the effective mass of a plate depends on its shape! 

• (2) For Lu S> c/c, there is dissipation due to a frequency dependent anisotropic shear viscosity, i.e. a type of 
'friction' in the vacuum. 

• (3) A second plate at a separation H modifies the mass renormalization by a function of kH, but does not 
change the dissipation for frequencies < (cfc)^ -I- {ttc/H)'^. 

• (4) For all frequencies higher than this first optical normal mode of the cavity, the mechanical response is 
infinite, implying that such modes can not be excited by any finite external force. This is intimately connected 
to the resonant particle creation reported in the literature [ P-Syiq ] . It is also an example of radiation from neutral 
bodies, a purely quantum effect. 

• (5) A phase angle 9 between two similarly corrugated plates results in Josephson-like effects: A static force 
proportional to sin(0), and an oscillating force for a uniform relative velocity. 

• (6) There is a (minute) correction to the velocity of capillary waves on the surface of mercury due to a small 
change in its surface tension. 

The rest of the paper is organized as follows. In Sec.^ we develop the formalism to calculate the effective action for 
time dependent small deformation fields. In Sec.|l| the behavior of the mechanical response functions are examined 
in the frequency-wavevector domain. We discuss the important special case of lateral oscillations of rough plates in 
Sec.0, and Sec.0 gives the conclusions of the paper, and suggests avenues for future explorations. In Appendix ^ 



the path integral method is applied to the standard Casimir effect due to the electromagnetic field, while some details 
of the calculations are included in Appendix ^ 
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II. PATH INTEGRAL FORMULATION OF THE DYNAMIC CASIMIR EFFECT 



In this section we introduce a quantization method for a field, subject to boundary conditions on surfaces which 
undergo dynamic shape fiuctuations. Although one can tackle the gauge invariant electromagnetic vector field by this 
formalism (see Appendix ^ , we use for simplicity the "scalar electrodynamics" model used most commonly in the 
literature [ [l0|j2q |. We consider a scalar field described by the classical action (c = 1) 



S = 



(3) 



where summation over fi ~ I, ■ ■ ■ , d in a d dimensional space-time is understood, and we have made a Wick rotation 
by introducing the imaginary time variable X'^ = it. We want to quantize the field subject to the constraint 
that it vanishes on a set of manifolds (objects) embedded in space-time. The manifolds are described by functions 
Xf^ — X^(yct); a D^-dimensional manifold embedded in d-dimensional space is parametrized by j/q = {Ua, ■ ■ ■ , Va")- 
Following Ref. the partition function of the system is written as 



(4) 



Z = ^ f fills {cf>(XM)) cxpL^S[A, 

= ^J 'Dcl^{X)flVi;^{y^) expl-^S[4>]+iY, J dy^ Va(j/a)</' (Xa(ya)) I , 

^ Q — 1 t <^ J 

where Zq is the partition function for the system with no boundaries, and ga is the determinant of the induced metric 
on the a'th manifold lEOl, 



9a, 



dX^ dX^ 



(5) 



required to make the integration measure over the manifolds reparametrization invariant. 

The integration over the field (f) can be performed, and the resulting expression for the partition function is 



where 



Z = J YlT>ija{ya) exp {-SillpaiVa)]} , 



Si[lp] = Ip'^Mlp = X! / '^yaV9^'^yfiV9l3'^o.{ya)G'^{Xa{ya) - X fj{y f)))^ p{y , 



a,f3 



and G'^{X — X') — {—d^d^)^^,. Integrating over the Gaussian fields {V'a} gives the final result of 



InZ 
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lndet{M(X„(y„))}. 



(6) 



(7) 



(8) 



We next focus on the specific example of two D (space-time) dimensional surfaces in d = Z) + 1 space-time with 
average separation iJ, and small deformations. The surfaces are thus parametrized by Xi{n) = (x, ft,i(x)) and 
X2(x) = (x, -I- /i2(x)), where /ii(x) and /i2(x) are the deformations. One can read off the matrix M from Eq.(|^) as 



M(x,y) 



VsTW G^i^ - y, /ii(x) - /ii(y)) Vg^^gAj) G^{x -y,H + h^i^) - h,{y)) 

V^V^ G^'i^ ~y,H + h2iy) - /ii(x)) G'^ix - y, h^i^) - h^iy)) 



(9) 



in which \/ ga (x) — ^1-1- {Vha)"^ where the gradient is with respect to the plane (space-time) coordinates (see 
Eq.(||)). Following Ref. |^, we perform a perturbative expansion of the matrix in terms of the deformation fields, as 



M(x,y) = Afo(x,y)-)-5A/(x,y). 



The matrix 



Mo(x,y) 



G'^(x-y,0) G'^(x-y,i7) 
G'*(x-y,i/) G'*(x-y,0) 



(10) 



(11) 



3 



describes two flat plates and is only a function of the differences of the coordinates. Hence it can be transformed 
using a Fourier basis into 



Mo(p,q) = 



G"^(p) G"^(p,i?) 



(27r)^<5^(p + q), 



where 



C^ip) J d^x G''{x,0) e'P'^, 
G"*(p, H) = J d^x G"'(x, H) e'P '', 



(12) 



(13) 



are the Fourier-transformed Green's functions. 

The part of the partition function that depends on the deformation fields can be written as In = — i Indetjl + 
AIq^SAI}. Expanding everything in powers of /i(x), and keeping only up to second order terms, we obtain (for details 
of the calculations see Appendix ^ 



lnZ„ = i J d^xd^y/^(x-y) [/ii(x)/ii(y) + /i2(x)/i2(y)] 
d^xd^yg(x-y) [/ii(x)/i2(y) + /^i(y)ft2(x)]. 



(14) 



The kernels in the above equation are defined as 

X(x) = afG(x,0)Fi(x) + Fi(x)i;^5(x) + F4(x)F6(x), 

g(x) = a2^(x,7?)i^4(x) +F2(X)2 +F3(X)2, 



where 



Fi(x 

F2(X 
J^3(X 
^^4(X 
^^5(X 
^6(X 



d^p G"^(p) 

d^G;H^ dG''{p,H) 
{2tt)D AA(p) dH 
d^p G'^{p,H) dG'^{p,H) 



(15) 



(16) 



dH 



{27r)D AA(p) 

_(Pp_G^HpJi) 
{2tt)D Af{p) ^ ' 

d^p G'^ip) fdG'^{p,H) 



ipx 



ip X 



(27r)^ Af{p) V 
d^p G'^{p,H) fdG'^{p,H) 



and 



(27r)^ Af{p) V '9^ 



AA(p)= [G'^(p)]'- [G'^(p,i/) 



ipx 



(17) 



We now focus on the case of interest in the Casimir problem which corresponds to c? = 4 and D ^ 3. The Green's 
functions can be calculated from the definition following Eq.(^, and using Eqs.(p^ as 



G^P) 

gHp,h) 

92g4(x,0) 
djG\x,H) 



1 

2^' 
1 

— e 
2p 



(18) 



-pH 



1 1 

2^ ^' 
1 - 3i/2 



27r2 {x^ + H^y 
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To obtain the behavior due to the time dependence of shape fluctuations of the surfaces, we rotate back to "real" 
time variables, changing Eq.jl^) to 

InZ,,^ - ^J dtdt'd'^xd'^y K{x~y,it~it') [hi{x,t)hi{y,t') + h2{x,t)h2{y,t')] (19) 

+ dtdt'd^xd^y Q{x ~ y,it~ it') [hi{x,t)h2{y,t') + hi{y,t')h2{x,t)\. 

One can diagonalize the above expression using the Fourier transformations 

K{x,t) = j e— *+^«- K{qM, (20) 

for the fields. The corresponding transformations for the kernels result in 

/C(g,w)= [ dtd^xe-"^'+"^-^K{x,it) (21) 



— oo 



dr / d X e 



and 

Q(g,w) = J dtd^xe-''^'+''^ ''Qix,it) (22) 

dr / d^x e~'^^+"i-''Q (^Vx^ + t^) 
= -i A^{q,uj), 

where we have performed a rotation of the integration contours from it to t. The symmetric appearance of the 
argument -s/a;^ + in Eqs.(^l]) and (|2^) is due to the underlying Lorentz invariance of the theory. Finally, the 
resulting expression for the effective action, defined via \nZ = iScs/ti, reads 

^ y / 7$ i^+il^^) {\hliq,^r + \h2iq,io)n (23) 
-A-{q,Uj) {hi{q,uj)h2{-q,-uj) + hi{-q, ~uj)h2{q,uj))} ■ 



III. ANALYTIC STRUCTURE OF THE RESPONSE KERNELS 



The kernels A± (g, lo) a re closely related to the mechanical response of the system. By substituting the expressions 
of Eq.(p^) in Eqs.(15-1^ and ( pl| , p^ ), after some manipulations, one obtains 



+ C30 



A±{q,uj) 



ds 



2^q^ -oj^/c^Hs/tt 



g±{s), toYUj<cq. 



(24) 



For LO > cq, the result is obtained by analytic continuation of the above. The need for such analytic continuation is a 
distinction between Minkowski and Euclidean spaces, and did not occur in Ref. [|. In Eq.(||), 



5±(s) 



1 ^ cosh s 
sinh^ s 



smh^ s 
cosh^ s 



± 2s 



sinhs 
cosh'^ s 



37rV4 



(s2 +7r2/4)2 



(25) 



are two functions of the dimensionless parameter s, and upper and lower cut-offs for the space-time are understood 
in case of divergences. It can be seen that the kernels are functions of the separation H , but depend on q and uj 
only through the combination = q^ — uP' jc^ . Rather than elaborating on the closed forms of the kernels, we shall 
describe their behavior in various regions of the parameter space. 
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The kernels can be calculated exactly for a single mirror. It is simpler to take the H ^ oo limit in the original 
Eqs.(^5Hl|), and Additional manipulations lead to A°^{q,u)) = 0, and 

A^{q,oj)=\ (26) 
[ ^^^{co^ - c^^ff^ brLO>cq, 

where sgn(a;) is the sign function. While the effective action and the corresponding response are real for > 0, 
they become purely imaginary for < 0. The analytic continuation to this regime is not unique. However, if we 
require causality for the response, we are led to the choice in which the imaginary part is an odd function, leading to 
the inclusion of the sign function in the above formula. An imaginary response function signals dissipation of energy 
psf , presumably by creation of photons as shown in the next section For the special case oi q = 0, Eq.(|2q) 

agrees with the results obtained previously for flat mirrors. (Note that Ref. considers the electromagnetic 
field rather than a scalar field, accounting for the discrepancy in the numerical prefactor. We have also done the same 
calculations in 1+1 dimensions and have reproduced the result of Refs. jlljjlj], who also use scalar fields, with the 
correct coefficient.) 

In the presence of a second plate (i.e. for finite H), the situation is more complicated. For > 0, we can examine 
the behavior of the integrand in Eq.(p^ and find that the kernels are real and finite in this regime. However, for 
< the situation is different. Applying contour integration, we can transform the integral in Eq.(|2j) to the form 

X = 27rj ^ residues - CI, (27) 

where CX is the contribution from the semicircle at infinity. The residues can be easily shown to be finite and real, 
and add to up to a finite and well behaved sum. Examining the contributions from the semicircle at infinity, we 
find terms with both real and imaginary parts, that generically behave as CT ^ exp[(_fi' — 2)L/ H]/[K{L/ H)^], with 
K = 2Q' H/tt, Q' ~ iQ. L is the radius of the semicircle, presumably corresponding to an upper cut-off in space-time. 
Now we can see that in the limit L — > oo, the situation drastically changes at i^T = 2. While CI vanishes for K < 2, 
it diverges with both real and imaginary parts for K > 2, where the inequality is understood in its strict sense. 

To summarize, we can subdivide the parameter space of the kernels into three different regions, as depicted in 
Fig. 1. In region I {Q^ > for any H), the kernels are finite and real, and hence there is no dissipation. In region Ila 
where — tt^/T?^ < < 0, the iJ- independent part of is imaginary, while the i?-dependent parts of both kernels 
are real and finite. (This is also the case at the boundary = —tt'^/H'^.) The dissipation in this regime is simply 
the sum of what would have been observed if the individual plates were decoupled, and unrelated to the separation 
H. By contrast, in region lib where < — vr^/i/^, both kernels diverge with infinite real and imaginary parts. This 
i?-depcndcnt divergence extends all the way to the negative axis, where it is switched off by a prefactor. 
Note that some care is necessary in the order of limits for (L, H) —f oo. 

Divergence of the kernels might be argued to signal a breakdown of the perturbation theory. However, one should 
note that the exponential divergence of the kernels is an essential singularity that can not be "cured" by standard 
renormalization techniques. It is in fact an artifact of the unphysical assumption of perfect reflectivity for the mirrors 
at all frequencies. Any natural frequency cut-off for the reflectivity of the mirror will round off this divergence p7[ , 
leading to a finite mechanical response. In other words, the imperfection in reflectivity at low frequencies provides 
a leakage mechanism for the cavity that stops the resonant energy build-up. Exponential creation of photons (with 
observation time) for perfectly reflecting boundaries has also been reported in the literature for the case of a one 
dimensional cavity ||l^ , using a completely different approach. 

IV. CORRUGATED PLATES 

We now concentrate on a concrete example, and examine the lateral vibrations of surfaces with fixed (time indepen- 
dent) roughness, such as two corrugated plates. We assume that the first plate undergoes a lateral motion described by 
r{t), while the second plate is stationary. The deformations of the plates are thus described by /ii(x, t) — hi{x~ r{t)) 
and /i2(x, t) = /i2(x). The lateral force exerted on the first plate is obtained from fi(t) = SScf[/6ri{t). Using Eq. (p3[) , 
and within linear response, this is given by 

Mu;)^X^A^)r,{cJ) + f^{cJ), (28) 

where the "mechanical response tensor" is 
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X^J{io) = hc J ^ q,q^ i^[A+iq,co) - A+{q,0)] (q)!' + (g, 0) (/^i (q)/i2 (-q) + /u (-q)/i2 (q)) | , (29) 
and there is a residual (static) force 

f^iu;) = 2^5{u:) J ^ iq.A^iq,0) (/^i (q)/i2 (-q) - /^l(-q)/i2(q)) ■ 



We can calculate the dissipation rate 

T/2 



(30) 



-T/2 

+ 00 

= ^ / 2^ ''^ ^^^■('^^ ni-i^)rjicj), 

— OO 

+ 00 

= / 2^ ^ IniXy (^) '"»(-'^)''i('^)> 

— oo 

in which we have used the fact that real (imaginary) part of a response function is an even (odd) function of frequency. 
Using the above formulas, we can predict various interesting effects for corrugated plates, that appear due to quantum 
fluctuations of vacuum. 



A. Single plate phenomena 



For a single corrugated plate, whose deformation is defined as /i(x) = dcosk-x, one can easily calculate the 
response tensor using the explicit formulas in Eq.(|2^). In the low-frequency limit, i.e. when lu ^ cfc, we can expand 
the result in powers of w. This gives Xij — SuiijUi'^ + 0{uj'^), where Sniij = AhcPk^kikj /{288tt^c), and can be regarded 
as corrections to the mass of the plate |Q. (Cut-off dependent mass corrections also appear, as in Ref. [^.) Note 
that these mass corrections are anisotropic, with 

Sm\\ Ahk^(f/{288TT^c), (32) 
Sm± = 0. 

Parallel and perpendicular components are defined with respect to k, and A denotes the area of the plates. The mass 
correction is proportional to h, and hence inherently very small. For example, if we have a macroscopic sample with 
d Ki X = 2iT/k « 1mm, density « 15gr/cm^, and thickness t Ki 1 mm, we find 8m/m ^ 10^'^"'. Even for deformations 
of a microscopic sample of atomic dimensions (close to the limits of the applicability of our continuum representations 
of the boundaries), 5m /ra can only be reduced to around 10"^". While the actual changes in mass are immeasurably 
small, the hope is that its anisotropy may be more accessible, say by comparing oscillation frequencies of a corrugated 
plate in two orthogonal directions. 

For oj 3> cfc, which is the high-frequency limit, the response function is imaginary, and we can define a frequency 
dependent effective shear viscosity by Xiji'-^) — —i^Viji^)- This viscosity is also anisotropic, with 

7?ll(w) = ?iAfc2d2,^V(7207rV), (33) 
11 d^) = 0. 

Note that in 1+3 dimensions the dissipation is proportional to the fifth time derivative of displacement, and there 
is no dissipation for a uniformly accelerating plate. However, a freely oscillating plate will undergo a damping of 
its motion. The characteristic decay time for a plate of mass M is t k, 2M/rj. For the macroscopic plate of the 
previous paragraph, vibrating at a frequency of w « 2cfc (in the lO^^Hz range), the decay time is enormous, r ^ lO^^s. 
However, since the decay time scales as the fifth power of the dimension, it can be reduced to 10~^^s, for plates of 
order of 10 atoms. However, the required frequencies in this case (in the lO^^Hz range) are very large. Also note 
that for the linearized forms to remain valid in this high frequency regime, we must require very small amplitudes, so 
that the typical velocities involved v ~ row, are smaller than the speed of light. These difficulties can be somewhat 
overcome by considering resonant dissipation in the presence of a second plate. 
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B. Double plate phenomena 



With two plates at an average distance H, the results are qualitatively the same for frequencies less than the natural 
resonance of the resulting cavity. There is a renormalization of mass in region I, and dissipation appears in region 
Ila, of Fig. 1. However, the mass renormalization at low frequencies {uj <^ ck) is now a function of both k and i7, 
with a crossover from the single plate behavior for kH '--^ 1. In the limit of kH <C 1, we obtain 

5m\\ = hABk^(f/48cH^, (34) 
Sm± — 0, 



with 

oo 

4 



B= / ds s 



2 







-0.452448. (35) 



Compared to the single plate, there is an enhancement by a factor of {kH)~^ in 5m\\. The effective dissipation in 
region Ila is simply the sum of those due to individual plates, and contains no H dependence. 

There are additional interesting phenomena resulting from resonances. Wc find that both real and imaginary parts 
of A±{q,Lo), diverge for up' l(? > + 7r^/if^, that corresponds to region lib of Fig. 1. In the example of corrugated 
plates, we replace q by fc to obtain a continuous spectrum of frequencies with diverging dissipation. (Note that the 
imaginary part of the response function gives the dissipation.) Related effects have been reported in the literature 
for 1+1 dimensions [|l5| - |lq| , but occuring at a discrete set of frequencies ujn — mrc/H with integer n > 2. These 
resonances occur when the frequency of the external perturbation matches the natural normal modes of the cavity, 
thus exciting quanta of such modes. In one space dimension, these modes are characterized by a discrete set of 
wavevectors that are integer multiples of -k / H . The restriction to n > 2 is a consequence of quantum electrodynamics 
being a 'free' theory (quadratic action): only two-photon states can be excited subject to conservation of energy. Thus 
the sum of the frequencies of the two photons should add up to the external frequency . 

In higher dimensions, the appropriate parameter is the combination up' jc? — cp . From the perspective of the 
excited photons, conservation of momentum requires that their two momenta add up to while energy conservation 
restricts the sum of their frequencies to w. The in-plane momentum g, introduces a continuous degree of freedom: the 
resonance condition can now be satisfied for a continuous spectrum, in analogy with optical resonators. In Ref. [ p7| , 
the lowest resonance frequency is found to be 27rc/iJ which seems to contradict our prediction. However, the absence 
of ui\ — t:c/H in 1-1-1 D is due to a vanishing prefactor |]l7| , which is also present in our calculations. However, in 
exploring the continuous frequency spectrum in higher dimensions, this single point is easily bypassed, and there is a 
divergence for all frequencies satisfying uP /P > + tP jYP' , where the inequality holds in its strict sense. 

Resonant dissipation has profound consequences for motion of plates. It implies that due to quantum fluctuations of 
vacuum, components of motion with frequencies in the range of divergences cannot be generated by any finite external 
force] The imaginary parts of the kernels are proportional to the total number of excited photons ||l^. Exciting 
these degrees of motion must be accompanied by the generation of an infinite number of photons; requiring an infinite 
amount of energy, and thus impossible. However, as pointed out in Ref. [^|, the divergence is rounded off by assuming 
finite reflectivity and transmissivity for the mirrors. Hence, in practice, t : 
the degree of ideality of the mirrors in the frequency region of interest. 



re restriction is softened and controlled by 



C. Josephson-like effects 



Consider two plates corrugated at the same wavelength, and separated at a distance H, whose deformations are 
described by /ii(x, t) = c?i cos[k • (x — r(i))] and ft.2(x, t) = c?2 cos[k • x]. For the case of uniform motion of the plate, 
namely, r{t) ~ v t + rg, we can find the full expression for the lateral force as 

m = ^ = ^A_(fc,0)kdid2sin[k.r(0]. (36) 
or[t) 2 

While in previous expressions we performed a expansion to linear order in r(i), in Eq.(p6[) the complete expression 
is calculated using the time translational invaraince of a uniformly moving object. For two stationary plates (v = 0) 
with a phase mismatch of a = k • rg, there is a (time independent) lateral force 
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Fdc = ——A^{k,0)kdid2sma, 



(37) 



which tends to keep the plates 180 degrees out of phase, i.e. mirror symmetric with respect to their mid-plane. This 
should be regarded as a correction to the static Casimir force due to the deformations of the plates. The dependence 
on the sine of the phase mismatch is reminiscent of the DC Josephson current in superconductor junctions, the force 
playing a role analogous to the current. (The amplitudes di and d2 are similar to the wave functions in SIS junctions.) 
Note that we could have equally well obtained this result using the constant term in Eq.(|30|). 

For uniform motion, there is an analog for the AC Josephson effect, with velocity (the variable conjugate to force) 
playing the role of voltage. Setting Tq = 0, we find 

hcA 

Fac = — A_(fc,0)kdid2sin[(k- v) t], (38) 

where A-{k,0) = g{kH), with g a (dimensionless) slowly varying function of its argument. Note that the force 

oscillates at a frequency w = k • v. Actually both effects are a consequence of the attractive nature of the Casimir 
force. It would be difficult to separate them from similar forces resulting from say, van der Waals attractions. One 
difference with (non-retarded) van der Waals force is in the power-law fall-off. Another potential difference with 
additive attractive interactions may be in the angular dependence of a static force. Such a difference in angular 
dependence due to the non-additive nature of thermal fluctuation- induced forces is discussed in Refs. pl|,p[ for two 
rods on a membrane. It is also likely to occur if we examine two rod-shaped deformations on the two plates in the 
Casimir geometry. 



D. Capillary waves 



In addition to solid corrugated plates, we can examine deformations of fluid surfaces. For example, consider the 
capillary waves on the surface of mercury, with a conducting plate placed at a separation H above the surface. The 
low frequency-wavevector expansion of the kernel due to quantum fluctuations in the intervening vacuum starts with 
quadratic forms and w^. These terms result in corrections to the (surface) mass density by 5p = hB/AScH^, and 
to the surface tension by Sa = hcB/48H^ {B is the numerical constant calculated in Eq.(|35|)). The latter correction 
is larger by a factor of (c/cg)^, and changes the velocity Cs, of capillary waves by Scs/c^ = ficB/96aH^, where tr is 
the bare surface tension of mercury. Taking H ^ 1mm and a ~ 500 dynes/cm, we find another very small correction 
of 5c,/cO ~ 10~i9. 



V. CONCLUSIONS AND OUTLOOK 



In conclusion, we have developed a path integral formulation for the study of quantum fluctuations in a cavity 
with dynamically deforming boundaries. As opposed to previous emphasis on spectra of emitted radiation, we focus 
on the mechanical response of the vacuum from which we extract a variety of interesting mechanical effects. Most 
of the predicted dynamic Casimir phenomena, while quite intriguing theoretically, appear to be beyond the reach of 
current experiment: the most promising candidates are the anisotropy in mass, and resonant dissipation. A major 
difficulty is the necessity of exciting deformation modes of macroscopic objects in the lO^^Hz or higher frequency 
range. However, as the static Casimir effect seems more amenable to current experiments Q, it would be interesting 
to test the non-additive nature of this force for plates with rod-shaped deformations. 

The path integral method is quite versatile: It is better suited to treating deformed boundaries than other techniques, 
such as that of Ref. |l^. Moreover, unlike other approaches it provides the dispersive part of the response, as 
well as its dissipative part. Furthermore, we can envision a number of future extensions based on this formalism, 
some of which are listed below: 

1. It is worthwhile to repeat the calculations for electromagnetic gauge fields. While we expect the above results 
to remain qualitatively unchanged, the vector nature of EM field, as well its gauge symmetry, may lead to novel 
additional phenomena. The transversality of photons in vacuum, as well differences between transverse electric 
(TE) and transverse magnetic (TM) modes, are two examples of the difference with scalar field, and may be 
important when considering the spectrum of emitted radiation (see below). 

2. There have been some efforts to calculate the spectrum of radiation emitted from a single fluctuating mirror, or 
by a cavity with fluctuating walls |0,E2[. Our formalism allows calculations in the more realistic case of 3+1 
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dimensional undulating walls. (Presumably, it is simpler to excite surface deformations of the boundaries at 
high frequencies, rather than rigidly oscillating the wall as a whole.) Radiation spectra can be obtained from 
two-point correlation functions which may be calculated by extending the methods of this paper. Hopefully, the 
angular distribution of the spectrum will provide signatures that can help in the experimental detection of such 
weak radiation. 

3. New issues can be probed by extending the calculation to the next order in the deformation fields. For example, 
new dissipative terms may arise at the second order, resulting in a 'viscosity' proportional to the velocity-squared 
as in Eq.(|). 

4. As in Ref. ||^, it is possible to examine other geometries such as fluctuating 1 and 2 dimensional manifolds, 
corresponding respectively to world-lines of a point particle and a wire in space-time. We can then use our 
techniques to calculate the dynamic Casimir forces involved for fluctuating wires, and mutual interactions 
between moving particles, wires, and surfaces. 

5. Dielectric super lattices are artificial media with nontrivial optical properties [Q. Because of the dielectric 
contrast in photonic crystals, they are excellent candidates for Casimir effects. In particular, dimensional 
arguments suggest that a photonic crystal with lattice spacing A, will receive a density enhancement (presumably 
anisotropic) of the order of ?i/(cA^). While this is at best many orders of magnitude smaller than a typical density, 
it is none-the-less interesting to calculate the associated corrections to the mass tensor. It is more promising 
to examine the coupling of the vibrations of such superlattices to their optimal properties. In particular, 
do mechanical deformations (phonons) modify the optical function of a photonic crystal, due to quantum 
fluctuations of the vacuum? 

6. An important generalization is to replace the assumption of perfectly metallic walls with boundary conditions 
that describe more realistic dielectric functions. 



ACKNOWLEDGMENTS 



RG acknowledges many helpful discussions with M.R.H. Khajehpour, B. Mashhoon, S. Randjbar-Daemi, and Y. 
Sobouti, and support from the Institute for Advanced Studies in Basic Sciences, Gava Zang, Zanjan, Iran. MK is 
supported by the NSF grant DMR-93-03667. 



APPENDIX A: PATH INTEGRAL METHOD FOR THE ELECTROMAGNETIC FIELD 

In this Appendix we introduce a path integral method for the standard static Casimir effect due to quantum 
fluctuations of the electromagnetic (EM) field. We start with the action 

So^-\J d'XV^ g^'^g'^^F.^F^p, (Al) 

where F^i, = d^A„ — d^A^,, and A^j, = {Ao,~A) is the vector potential. We have introduced the constant (space-time 
independent) metric as a source, to enable calculation of the integral over the stress-energy tensor |Q, as 

= 1 fd^XT,,{X), (A2) 
jji^f ^ J 

and 7]'^" = diag(— 1, 1, 1, 1) is the flat space-time metric. 

We then take two perfectly conducting parallel plates separated a distance H along the z-direction. The boundary 
condition subject to which we have to quantize the theory is that the transverse electric field should vanish on the 
plates 1^. We should also treat the gauge freedom of the action by adopting a proper gauge fixing procedure. One 
can show that a possible proper quantization is achieved by considering the following expression for the generating 
function p5[ 

Zig) = lim — / VA{X) TT TT S {doAiXM) - d^AoiX^ivo.))) e^^'^I^l, (A3) 
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where 



EE So[A] - Ay d^xV^ig'^'^dM' 



(A4) 



1 



Now one can easily calculate the average value for energy using 

d\nZ{g) 



2ih 



dg' 



,00 



d^X (TooiX)) 
ATS, 



(A5) 



where A is the area of the plates, T is the overall time interval, and £ is the Casimir energy. We now represent the 
delta-functions using Lagrange multiplier fields, as 



Z{g)= lini^ / VA{X) V^,^^{y^) exp <j ^Sp[A]+iJ2 I dy^ V'»,a(ya) (5oA:(^a(ya)) - 9»v4o(^a(ya))) 



Performing the integration over the gauge field, yields 



1 

Z{g) = lim 1^ / TT Vlpi^aiVa) exp{-S'[V'j,a(2/a)]}, 



(A6) 



(A7) 



in which 



+do'ipi,a{ya)Gij{Xa{ya) - Xp{y0))do1pj^a{ya) 
~ 2doll;i^a{ya)Goi{Xa{ya) - Xf3{yf3))dj-ipj^a{ya)] 



(A8) 



and the matrix Green's functions satisfy 



G,.x{X ^ X') = Si^^S^X - X'). 



(A9) 



The geometry of two parallel flat plates at a separation H is implem ente d using the embeddings Xi(x, t) = {t, x, 0), 
and X2(x, t) = {t, X, H). Then we can read off the matrix M from Eq.( A8). In Fourier space, it is diagonal and reads 



M(p, q) = (27r)^(5^(p + q) {p^PjMao + pIm,^ + poPiMoj + poPjMo,) , 



in which 



and 



Gm.(p,o) g^4p,h) 

G^,{p,H) G^,(p,0) 



G^Up,H) = J dtd^^G^Jt,^,H) e^P— ^Po*. 



Having found the matrix M, we can calculate the generating function from 

In Z{g) = -ilndet{Af}. 



(AlO) 



(All) 



(A12) 



(A13) 



To this end, we only need to determine the matrix Green's functions from Eq.( ]AS| ), which in momentum space are 
given by 
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9u\ 



g'^'PcPp 



(A14) 



where Pi, is the four- momentum. At this point we can set [3 = 0, which implements the (Lorentz) gauge choice 
without any difficulty. 

To simplify our calculations, we take the specific form g'^" = diag((7°", 1, 1, 1) for the metric. This is sufficient for 
our purposes, as we are only need variations with respect to g^" to obtain the average energy. Inserting the matrix 
Green's functions from Eq.(Al4), in Eqs.( AlO- A12| ), we obtain after some manipulations. 



Mip,q) = i2nrS'{p + q) 



PiPj 



/(5°°,p,0) /(.gO",p,i/) 
I{g°",P,H) /(.g"o,p,0) 



where 



I{9'",P,H) 



We can next calculate the determinant to get 

d^p 



InZ(g) 



-AT 



(27r)3 



Finally, we use Eq.(A5), and obtain 



exp[-g(ff>g + p^)V^] 
2{g00pl + p2y/2 ■ 

ln(l-exp [-2i/((?>2_^p2)i/2 



TT^ he 

720^' 



which is the well known result for the Casimir energy |^ . 



(A15) 



(A16) 



(A17) 



(A18) 



APPENDIX B: PERTURBATION AROUND FLAT PLATES 



In this Appendix we present the details of the calculations that lead to Eq. (|T^) . We follow the steps and notations 
in the Appendix of Ref. j^. We also generalize this work by including deformations of the second plate, and by using 
an invariant measure for integrations over the manifolds. We also note that the final result of Eq.(14), corrects some 
mistakes that appeared in Ref. To find 6M we expand the square roots, as well as the Green's functions in Eq.(||) 
up to second order in the deformation 



Bi(p,q) A(p,q) 
yl(q,p) B2(p,q) 



G'^(x - y, /ii,2(x) - h^^2{y)) = G'^(x - y, 0) + -d^G^i^ - y, 0)[/ii,2(x) - /.l,2(y)]^ 
G'^(x -y,H + /i2(x) - h,{y)) = G'^(x -y,H)+ d^G^i^ - y, i?)(/i2(x) - h,{y)) 

+ ld!G'i^-y,H)[h,i^)-h,iy)f. 

After Fourier transformation one obtains (see Eq.([Tc|)) 

JM(p,q) 
where 

A(p,q) = J d^xd^y e^P-^+^^ y I^^G\^-y,H) [{Wh,{^))' + {Wh^iy))'] 

+d,G''{^-y,H){h2{y)~h,{^)) 
+ ia2G'^(x-y,H)[/i2(y)-/Ji(x) 

Bi,2(p,q) = J d^xd^y e'P-+''i '' |iG''(x-y,0) [{Vh,,2{^))' + i^hM)'] 

-fi^f G'^Cx - y, 0)[/ii,2(y) - /ii,2(x) 



(Bl) 
(B2) 



(B3) 



(B4) 
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The matrix Mq is diagonal in Fourier space and can be easily inverted as 

1 



^V'(p,q) 



Af{p) 



(27r)^<5^(p + q). 



(B5) 



We then calculate InZh = -5 trln(l + Mq ^5M), using the expansion of the logarithm up to the second order. We 
also choose / d^x /i(x) = without loss of generality, leading to 



InZf,, = -- 



d^^d^y{^c'P-^^-''4G'^(x-y,0)((V/.r(x))^ + (V;^.(x))2) (B6) 
+ ia,2G'^(x - y, 0) ((/ii(y) - h.i^))' + {h^{y) - h^{^)f) 
-^W^ e^^-^-'') [G'^(x-y,i7) ((V/.i(x))^ + {Vh,{^)f) 

1 [ d^p d^q (G'^{p)G'^{ci)+G'\p,H)G'^{ci,Hy 



X d^xd^v e^'i-^+^P-^ d.G'^iy - x, i7)/i2(x) - j d^xd^v e^P'^+'i'^ 9.G"^(v - x, H)hi{^) 
X y d^yd^ue-"i y-'P "5,G'^(u-y,i7)/i2(y)-y" d^yd^u e-'P ™" a,G'^(u - y, ii-)/ii(y)| . 

Integration over u — y and v — x can now be performed followed by integrations over p and q, resulting in Eq. (^4|) . 
Note that we have eliminated some uninteresting terms that were proportional to J d^p and J d^p p^. Such integrals 
are usually removed in the framework of dimensional regularization ]3^ , and don't play any important role. They 
appear since the "constant" infinite energy that is usually neglected when quantizing a field, now depends on the 
configuration of the surfaces. When we dynamically deform the boundaries, this contribution changes. These terms 
give rise to cut-off corrections to mass and surface tension as discussed in Ref. [F9I 
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FIG. 1. Different regions of the (q, tij) plane. 
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